Critical dynamics of ballistic and Brownian particles in a heterogeneous environment 
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The dynamic properties of a classical tracer particle in a random, disordered medium are investigated close 
to the localization transition. For Lorentz models obeying Newtonian and diffusive motion at the microscale, 
we have performed large-scale computer simulations, demonstrating that universality holds at long times in the 
immediate vicinity of the transition. The scaling function describing the crossover from anomalous transport to 
diffusive motion is found to vary extremely slowly and spans at least 5 decades in time. To extract the scaling 
function, one has to allow for the leading universal corrections to scaling. Our findings suggest that apparent 
power laws with varying exponents generically occur and dominate experimentally accessible time windows 
as soon as the heterogeneities cover a decade in length scale. We extract the divergent length scales, quantify 
the spatial heterogeneities in terms of the non-Gaussian parameter, and corroborate our results by a thorough 
finite-size analysis. 
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I. INTRODUCTION 

Heterogeneous materials abound in synthetic products and 
in nature; they are composed of domains of different ma- 
terials or phases, with characteristic dimensions covering a 
wide range of length scales. A physical understanding of 
their macroscopic properties, such as mechanical elasticity, 
electrical conductivity, particle transport, or fluid permeabil- 
ity has far reaching consequences for applications in mate- 
rial science, nano-chemistry, oil recovery, and even biology. 
Examples include anomalous transport of tracers in porous 
soil columns slow diffusion of sodium ions in sodium 
silicates jSJjl, and transport in colloidal gels close to gela- 
tion |0, HlaS [Ml In biology, the dense packing of differ- 
ently sized proteins, lipids, and sugars in the cell cytoplasm 
is summarized as macromolecular crowding J^, [T(| Ell- It 
leads to a suppression of diffusion with increasing molecu- 
lar weight IU2I| and to anomalous diffusion fl3ll . observed in 
eukaryotes HI d [M El B 

and bacteria 1191 1. 

Transport of tagged ions, macromolecules, or nanoparti- 
cles in such heterogeneous environments is strongly hindered, 
since the presence of a variety of components reduces the 
accessible volume to a small fraction of three-dimensional 
space. In computer simulations, one observes a drastic sup- 
pression of the diffusion coefficient by up to several orders of 
magnitude upon decreasing the porosity, i. e., the fraction of 
the non-excluded volume. There, the medium is often mod- 
eled by randomly placed obstacles lf20L I2H l22tl . but recent 
studies investigate also realistically reconstructed media, e. g., 
Vycor glass JH] and North Sea chalk JM). 

In the above examples, three major transport phenomena 
are observed: normal diffusion, immobilization or localiza- 
tion, and anomalous transport. We will demonstrate in the 
following that all three aspects may be unified into the concept 
of transport in a disordered, heterogeneous medium with a 
percolation transition; such a transition entails a critical point 
with a divergent correlation length. The asymptotic behav- 
ior of this length scale in the critical region, together with the 
intrinsic properties at criticality, is encoded in the renormal- 



ization group flow; therefrom, all macroscopic observables 
(such as the diffusion coefficient) can be inferred in principle. 
This leads to the postulate of universality: systems sharing 
the same critical point exhibit the same universal scaling laws 
in the critical regime. Consequently, one expects a generic 
mechanism for slow, anomalous transport in a heterogeneous 
environment. 

In a recent study on the Lorentz model, i. e., for a ballistic 
tracer in a porous medium, we have shown that a continuum 
percolation transition of the accessible volume is responsible 
for the suppression of the diffusion coefficient ll25ll . We have 
successfully applied the theory of critical dynamic scaling to 
explain the dynamics over many decades of time and length 
scales, covering a large range of porosities. This analysis will 
be extended and further substantiated here. In addition, we 
present simulation results for Brownian particles in a porous 
medium, modeling macroscopic particles, e. g., proteins, that 
experience thermal fluctuations from a solvent. Universality 
implies that the long-time behavior of dynamic observables 
close to the transition is independent of the microscopic de- 
tails. It predicts that the critical exponents of the anomalous 
diffusion and the scaling of the diffusion coefficient are the 
same and moreover that the scaling functions coincide. The 
equivalence of ballistic and Brownian particles is well known 
for molecular and colloidal glasses [26], and it is expected to 
hold generally for slow dynamics. By a direct comparison of 
ballistic and Brownian particles, we test this hypothesis and 
give an estimate of the asymptotic scaling regime. The lat- 
ter is essential for the interpretation of experiments, since it 
allows to assess the applicability of the asymptotic laws to a 
specific measurement window. 

Let us summarize the different transport phenomena briefly. 
Diffusion is well known and understood since more than a 
century now ll27il28ll . It is connected with Brownian motion, 
observed for a large particle kicked around by the surround- 
ing gas or solvent molecules. The corresponding mean-square 
displacement grows linearly in time, 



Sr 2 (t) := (AR(t) 2 ) = 2dDt, 



(1) 
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where AR(i) denotes the displacement of the particle after 
a sufficiently long time lag t, d is the space dimension, and 
. . . ) an appropriate average. The phenomenon of localiza- 
tion is reflected in a plateau in the mean-square displacement, 
defining the localization length £, 

8r 2 (t)= const = t 2 . (2) 

Particles may be immobilized due to chemical binding, or they 
get trapped in cages formed by the surroundings. 

Anomalous transport is certainly the most fascinating phe- 
nomenon among the three. Contrary to normal diffusion, it 
is not simply a consequence of the central limit theorem; in 
particular, it requires the non-trivial presence of either a broad 
distribution of time scales or long-ranged correlations l29ll . 
The mean-square displacement exhibits a power law growth 
with a fractal exponent, 

Sr 2 (t) ~t 2/d "; (3) 

here we shall discuss subdiffusion, e? w > 2. In most cases, it 
originates in one of two mechanisms: binding to finite traps 
with a broad distribution of binding rates, or confined motion 
in a spatially non-uniform, heterogeneous medium. Both sce- 
narios may be easily distinguished experimentally, since bind- 
ing rates obey an Arrhenius law in general, and thus, the dy- 
namic coefficients are very sensitive to temperature changes. 
Steric interaction in contrast is insensitive to temperature and 
often conveniently modeled by hard potentials. 

The article is organized as follows: In the first part, we in- 
troduce the Lorentz model (Section |Hj, outline its relation to 
continuum percolation, and review some theoretical results 
relevant for the subsequent analysis (Section [TITb . In Sec- 
tion llVI we set forth the applied simulation techniques, while 
Section|V]is devoted to the scaling behavior of diffusion coef- 
ficients, length scales, the mean-square displacement, and the 
non-Gaussian parameter, complemented by a finite-size scal- 
ing analysis. 

H. THE LORENTZ MODEL 

The Lorentz model is a generic model covering all of the 
above aspects of transport in a heterogeneous environment. It 
has attracted the attention of researchers in statistical physics 
over more than a century by now and was applied to various 
contexts. Lorentz introduced the model in 1905 as a micro- 
scopic justification of the Drude conductivity of a metal ll3"oll : 
a non-interacting electron gas is scattered off randomly dis- 
tributed ions. Equivalently, one may think of a single ballistic 
and point-like tracer particle with velocity v = |v|, explor- 
ing the void space between randomly, uniformly, and inde- 
pendently placed scatterers of number density n. A hard-core 
interaction with the scatterers is usually employed, yielding 
impenetrable obstacles of radius a. This interaction poten- 
tial is equivalent to a tracer particle and obstacles sharing the 
same radius a/2. The kinetic energy of the particle is con- 
served, and the unit time and length scales are set by a and 
to := v~ x a, respectively. The only control parameter of the 



model is the dimensionless obstacle density n* := na d ; it is 
directly linked with the porosity of the medium, 

ip = exp(— 47rn*/3). (4) 

Statistical averages include different initial positions (re- 
stricted to the void space) and velocities (with fixed magni- 
tude) and different realizations of the disorder; in particular, 
the tracer particle is not restricted to the percolating cluster. 

Although the model is a deterministic system, the initial ve- 
locity of the particle is quickly randomized by subsequent col- 
lisions with the obstacles, yielding diffusive motion on large 
time scales. Such a stochastic trajectory may be thought of as 
a hopping process with a general distribution of waiting times 
p T (r) and displacements p a (a) l29ll . Provided that the first 
and second moments exist, uncorrelated steps generate diffu- 
sive motion on large time scales with diffusion coefficient 

(a 2 ) v 2 (t 2 ) 
D = -> — '— = — -> — '-■ (5) 

2d{r) 2d(r) ' ' 

the second form refers to a ballistic particle with fixed ki- 
netic energy. In the Lorentz model, the collisions constitute 
a Poisson process with average rate r c _1 , thus the waiting 
time distribution has first and second moments /r) = r c and 
(r 2 ) = 2r 2 , respectively. The mean collision rate follows 
from the condition to encounter no obstacles in a corridor 
of volume it a 2 (vr c ) = 1/ninrf = 3 dimensions, and thus 
t" 1 = irn*v/a. Since the differential scattering cross-section 
for a sphere is isotropic, it corresponds to the transport cross- 
section describing the transfer of momentum; subsequent col- 
lisions are uncorrelated in the dilute limit, n* — > 0. Under 
these conditions, Eq. (0 yields the correct diffusion coeffi- 
cient H, 

D a = —: for 7i*^0, (6) 
iirn* 

recovering the leading order in an asymptotic low-density ex- 
pansion Enm. 

At finite densities, spatial correlations between obstacles 
induce persistent anti-correlations of the velocity, which are 
reflected in a negative tail of the velocity auto-correlation 

function, ip(t) = (v(i) • v(0)) /v 2 t -d/2-i for f 

00 ll33l, l34ll . Since the integral over ip(t) is related to the dif- 
fusion coefficient, such a tail reduces D. The long-time tail 
can be thought of as a consequence of backscattering events 
preferring the return of the particle to its origin, emphasizing 
the importance of the topology of the obstacle matrix. The 
exponent of the tail is universal, i. e., it does not depend on 
the density of the scatterers. But as the density is increased, 
a pre-asymptotic, negative tail emerges, which suppresses the 
diffusion coefficient additionally I3~5t[36ll3~7ll . At a critical ob- 
stacle density n*, the pre-asymptotic tail persists for all times, 
yielding exactly zero diffusivity. 

At higher densities, diffusion is absent too: the tagged par- 
ticle is trapped, and the mean-square displacement saturates. 
The transition from diffusion to localization has the signatures 
of a continuous phase transition; it exhibits power-law diver- 
gences of physical quantities with universal critical exponents. 
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In particular, the diffusion coefficient vanishes with a power 
law upon approaching the critical density, 

D ~ |n* -n*| /1 . (7) 

The localization transition was predicted by a mode-coupling 
approa ch l35l |36l l38ll and by a mapping to continuum perco- 
lation 11391 14011 . The direct link between the Lorentz model 
and continuum percolation, however, was established only re- 
cently B25I1 : it is substantiated further by the data presented in 
Section [V] Since the interpretation of the data depends on an 
understanding of continuum percolation theory and some of 
its subtleties, the next section summarizes and motivates the 
central results of the theory. 



III. CONTINUUM PERCOLATION 

One particular interesting class of percolating systems is 
that of continuum percolation ||4lll. which was later also 
termed "Swiss cheese" model 14211 . Isolating discs or spheres 
are distributed at random within a conducting material, they 
may overlap, and the remaining material between them forms 
the percolation clusters — with respect to the medium, this 
model is identical to the Lorentz model. In the simplest vari- 
ant, the spheres are monodisperse, i. e., have equal radii, but 
binary distributions were considered as well IH31 Hp] , The 
percolation threshold was first determined by discretizing the 
void space and using lattice methods ll4lll . Much more effi- 
ciently, the void space is well represented by a random net- 
work constructed from a Voronoi tessellation l45tl : then, the 
percolation thresholds are calculated for that network. As ex- 
pected from the universality hypothesis, the critical exponents 
for static properties, e. g., f3 for the strength of the percolating 
cluster Poo and v for the correlation length £, agree with the 
findings for various lattices of the same dimension l46ll . 

The more suprising were theoretical predictions, based on 
the "nodes-links-blobs" picture, that the exponents of trans- 
port properties, e.g., conductivity and elasticity, are consid- 
erably larger than their universal values on lattices j40l l42tl 
(with the exception of the conductivity exponent in d = 2). 
These results were obtained by assigning physical proper- 
ties to the bonds of the network, e. g., in order to investi- 
gate the conductivity, the bonds were considered as resis- 
tors. Since the macroscopic conductivity only depends on the 
probability distribution of the bond conductances, the latter 
may be distributed randomly obeying a given distribution. A 
renormalization group analysis of such random resistor net- 
works II47[ l48l |49j |50J attributes the non-universal character 
of the critical exponents to a singularity in the distribution of 
conductances. In the following, we will sketch the derivation 
of such a singular distribution for the Lorentz model, and we 
will summarize the line of argument for the modification of 
the tracer dynamics. 



A. Mapping to random resistor networks 

The mapping from the Swiss cheese model to a geometric 
Voronoi network to a random resistor network is indeed intu- 
itive, but not at all rigorous. Further, the ballistic dynamics of 
the tracer particle in the Lorentz model is quite different from 
the diffusive dynamics of a random walker on a percolation 
lattice. Only by virtue of the universality hypothesis, one can 
expect that both models share the same critical behavior. 

For the dynamic properties of the three-dimensional Lo- 
rentz model, it is essential that narrow gaps between obstacles 
are abundant, i.e., the gap distribution exhibits a singularity. 
The mapping from the Lorentz model to a random resistor 
network includes the following steps (see Fig. [TJ: 

(i) Construct a Voronoi tessellation around the obstacle 
centers, defining vertices that are connected by edges. 
Each vertex either defines a "chamber" between the ob- 
stacles or is outside the void space. Each edge or bond 
defines a "gap" between three obstacles at a time. (The 
plane through the obstacle centers separates adjacent 
chambers.) If a gap is blocked by these obstacles, re- 
move the related bond. Then, any path in the void space 
can be transformed continuously int o a p ath along the 
bonds without crossing the obstacles IH5I1 . 

(ii) Choose a chamber and one of its gaps. A transition rate 
W can be assigned to this gap, given as the ratio of the 
phase-space volume available for leaving through the 
gap to the total volume V of the chamber. The numera- 
tor is proportional to the cross-section area of the gap A, 
and W oc A/V lf40h . Slow transport is connected with 
the distribution p(W) for small transition rates. Hence, 
the singular behavior of p(W) for W — > is deter- 
mined by the distribution of small cross-sections A or 
narrow gaps. In d dimensions, one calculates lf40h 

p(W) ~ W~ a , where a = (8) 

d — 1 

(iii) Drop the obstacles, keep only vertices and bonds in the 
void space, and interpret the transition rates through the 
gaps as conductances along the bonds. This defines the 
random resistor network with a power-law distribution 
of weak conductances. 

Such random resistor networks have been investigated exten- 
sively by means of Monte Carlo simulations IdH l52ll and 
renormalization group techniques ll48ll49Tl . providing reliable 
numeric and analytic results for the critical behavior f53ll . 

B. Conductivity 

The central property of random resistor networks is the 
(macroscopic) conductivity E. As the infinite cluster becomes 
increasingly thinner upon approaching the threshold, it van- 
ishes with a power law, E ~ |e| M , anticipating the same ex- 
ponent fj, as for the diffusion coefficient due to the celebrated 
Einstein relation D ~ E l53ll . 



FIG. 1: (Color online) Mapping of the Lorentz model to a random resistor network: (i) construct the Voronoi tessellation around the obstacle 
centers, (ii) remove blocked bonds and assign a conductivity to each bond according to the gap width, and (iii) finally drop the obstacles. 
(Cartoons are in two dimensions.) 



The percolating network looks homogeneous at length 
scales larger than a correlation length £. This fact is pictorially 
reflected in the "links, nodes, and blobs" model by Skal and 
Shklovskii JH] and de Gennes JH); for reviews see JUIH. 
It simplifies the infinite cluster of the percolating network in 
the following way (Fig. |2): the nodes are defined as vertices 
which can not be isolated from the cluster by cutting any two 
(not necessarily adjacent) bonds. The "nodes" are supposed 
to be distributed homogeneously with an average spacing of 
the correlation length £. Most vertices of the network belong 
to cul-de-sacs ("dangling ends"), i. e., they can be separated 
by cutting a single bond somewhere in the network. Nodes 
are connected by "chains", i. e., objects that can be isolated by 
cutting two bonds adjacent to a node. A fraction of the bonds 
of a chain carries the whole current ("links" or "red bonds"), 
the remaining bonds of a chain form "blobs" which connect 
two links at a time. 

The chains have an average resistance 1Z, which exhibits a 
power-law divergence with exponent ( upon approaching the 
percolation threshold, 

ft~H~ C - (9) 

Here we have introduced the separation parameter e; in case 
of the Lorentz model, we use e := (n* — n* c )/n* c . Apply- 
ing an electric field E, the potential drop between two nodes 
is of order yielding a chain current / = £E/1Z. Then, 
the current density in the network is given by j = 
Plugging in the definition of the conductivity, j = TjE, yields 
E ~ £ 2 ~ d /TZ, and a hyperscaling relation follows, 

li=(d-2)v + C (10) 

The remaining task is to identify the exponent £ for such 
random resistor networks that are relevant for the Lorentz 
model, i. e., that exhibit a power-law distribution of the con- 
ductances, p(W) ~ W~ a with < a < 1. We will follow 
the argumentation of Straley l47ll . The distribution of chain 
conductances pamn(W) can be considered as the renormal- 
ized distribution of bond conductances p(W) in the sense of 
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Dangling Ends 

FIG. 2: (Color online) Cartoon of the nodes-links-blobs model: 
nodes, being distributed homogeneously with average spacing £, are 
connected by chains. A chain is made of a series of links and blobs; 
links are defined as bonds carrying the whole current of a chain. Most 
bonds of the network are dangling sites, i. e., they carry zero current. 



the real space renormalization group. The renormalization 
flow shifts the whole distribution; usually, the peak of the 
renormalized distribution pchain (W) is determined by the peak 
of the microscopic distribution p(W). The low-conductivity 
tail, however, is shifted as well. Two different scenarios arise 
dependent on the value of a: either the peak absorbs the 
tail, yielding the universal form of the renormalized distri- 
bution, and the chain resistance exhibits the universal expo- 
nent £ unlv . In the second case, the tail pulls weight out of the 
peak, leading to a final distribution with the same small-M^ 
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tail, p c hmn{W) ~ W~ a . Then, the chain resistance is domi- 
nated by the small conductances. The total resistance of par- 
allel bonds, i. e., the blobs, may be neglected compared to the 
resistance of the red bonds, and p c hain(W) obeys for small W 



Pchain (W) ~ N led p(W), 



(11) 



where iV rei j is the number of red bonds within a chain. Dimen- 
sional analysis suggests the scaling form 



(12) 



with the typical chain conductance Wo- From Eq. dTTb one 
infers N lel j ~ Wq , and the average chain resistance TZ 

I 1 1' 1 is identified as Wo- For the divergence of the num- 
ber of red bonds, we finally employ a result by Coniglio HH], 



Collecting results, 



K~ W, 



o 



N. 



l/(a-l) 



(13) 



(14) 



yields £ = (1 — a) l . It has been argued that the crossover 
between both scenarios occurs such that l47ll59tl 



c 



(15) 



Without resorting to the assumptions of the "links, nodes, 
and blobs" model, Stenull and Janssen |50ll have proven di- 
rectly by means of an expansion of the renormalization group 
equations in e = 6 — dto arbitrary order the equivalent relation 



p = max [(d — 2)v + (1 — a) 1 , p 



Ian 



(16) 



where /j, lat is the universal exponent for lattice percolation. 

From simulations, it is known that 1 < £ unlv < 1.3 for 
d > 2 dimensions lf57tl . Using the above value of a for the 
Lorentz model, Eq. ©, it follows that £ takes its universal 
value on lattices only for d ~ 2 dimensions, otherwise £ = 
d-1. 



where the factor allows for the zero diffusion coeffi- 



cient of particles on finite clusters. Recalling D 



and 



£ ~ |e| v , the conductivity exponent is related to the walk 
dimension by 



p = (d w - T)v + (3. 



(19) 



In the above calculation, the diffusion coefficients of parti- 
cles from different clusters were averaged. Instead of taking 
the cluster average in the end, one may average already the 
mean-square displacement. Then, the relevant dynamic length 
scale is the root-mean-square cluster radius |25,[60ll. 



(20) 



and anomalous transport is characterized by the dynamic ex- 
ponent z different from d w l53ll . 



5r 2 (t)~t 2/z for £<< x . 



(21) 



The crossover time scale t x , being the same for all clusters, is 
not affected by the cluster average. It holds l z ~ £ d " ~ t x , 
and therewith D ~ £ 2 /i x ~ £ 2 ~ z , which yields a scaling 
relation connecting p and z, 



2v - (3 + p 



(22) 



In particular, this relation implies a connection of the expo- 
nents z and C, via Eq. ( tlO) . and z is expected to deviate from 
its universal value on lattices for d > 3. Moreover, z can be 
calculated from the geometric exponents v and 3 which are 
believed to equal their universal lattice values |46ll; in three 
dimensions, we use v = 0.88 and 8 = 0.41 15711 throughout 
this work, evaluating Eq. (l22l to z = 6.25. 



IV. SIMULATION DETAILS 



C. Transport on percolation networks 

Going back to the Lorentz model in the series of mappings, 
the transport of a particle in the percolating void space may be 
understood as a random hopping process between the nodes of 
the network with lattice constant £. Hence, the motion of such 
a particle is diffusive at scales much larger than the correlation 
length £. Note that £ also characterizes the linear dimension of 
the largest finite clusters. Along the fractal chains of the net- 
work, however, the dynamics is slow and exhibits anomalous 
transport characterized by the walk dimension d w > 2, 

< AR W 2 >chai„^ 2M "- (17) 

The time scale associated with £ obeys t x ~ £ dw , and diffusion 
is observed for long times, t 3> t x , with a coefficient accord- 
ing to Eq. ©. If the particle starts at any cluster, the diffusion 
coefficient is given by 

D=p 00 e/t*~\sfe- d ", (is) 



Molecular Dynamics simulations allow for a direct numer- 
ical analysis of the dynamic properties of the Lorentz model 
without resorting to random resistor networks. Thus, a quanti- 
tative description over the full density range becomes accessi- 
ble. Ballistic trajectories are produced by means of a standard 
simulation algorithm already employed by Bruin 16111 . It is 
combined with a method for calculating correlation functions 
online, optimized for exponentially large time scales. For 
Brownian particles, we have extended this simulation algo- 
rithm to include stochastic forces similar to recently discussed 
ideas JH. 

A. Ballistic particles 

We use an event-oriented simulation algorithm, since the 
tracer particle and the obstacles interact via a hard-core po- 
tential. The algorithm propagates the particle freely from col- 
lision to collision, in each case calculating the precise point in 
time of the next interaction with an obstacle. If the particle is 
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located at r with velocity v, it possibly hits a sing le obstacle 
in the coordinate origin after the time interval 16911 



toon = ~K~ \\fV-v 2 (r*-<j*) 



(23) 



where b = v-r. If b > 0, the particle departs from the obstacle 
and no collision will take place. If the radicand becomes neg- 
ative, the particle misses the obstacle. The particle is scattered 
specularly, its post-collisional velocity is 



v' = v — 2(v • &) <T, 



(24) 



where <r = (vi co n + r)/cr specifies the surface normal at the 
collision point. 

Therewith, the implemented algorithm is straightforward: 
In a preparation step, the obstacle positions and the initial 
phase space coordinates of the particle are chosen randomly 
from a uniform distribution, under the constraint that the par- 
ticle starts in the void space and |v| = v. To reduce the 
number of collision tests, the simulation box is divided into 
small cubic cells, on average containing one or two obstacles; 
only obstacles from the cell where the particle is located and 
from neighboring cells are considered. Then, repeated colli- 
sion tests and propagation of the particle to the next collision 
point yield the trajectory. 

Since a single trajectory can consist of several billion colli- 
sions, one easily runs out of computer memory if one would 
store the complete trajectory. An efficient blocking scheme 
("order-n algorithm", see [63]) takes care of this issue by ar- 
ranging the trajectory on a logarithmic time grid and simul- 
taneously calculating various correlation functions C(t;to). 
The algorithm already averages over different time origins to 
("moving time average"), which are, however, not necessarily 
uncorrected, especially for short time intervals t — to. Hence, 
we will not infer any estimate of the statistical error from this 
averaging procedure. Rather, a set of N t trajectories with dif- 
ferent initial positions for each of N r different realizations of 
the obstacle disorder is simulated, and the statistical error is 
estimated from Nt x N r independent measurements. At each 
density, we have simulated at least N r = 20 realizations of 
the disorder. At intermediate densities, the total number of 
trajectories has been chosen N t x N r > 100. This value has 
been increased up to 600 at very high densities, where the 
phase space is highly decomposed into small, disconnected 
parts. The longest trajectories span about 10 10 collisions, the 
demand on CPU time for such a trajectory was about 15 hours 
on a single AMD Opteron 248 processor core. 

All numerical results presented in this work refer to fixed 
dimensionality, d = 3. The simulation box has periodic 
boundaries, and its linear size was chosen as L = 200cr. A 
detailed finite-size analysis is presented in Section lVDl 



B. Brownian particles 

When the pores between the obstacles are filled with some 
solvent, the tracer particle performs Brownian motion on mi- 
croscopic time and length scales. By another Einstein relation, 




FIG. 3: (Color online) Typical particle trajectory in two dimensions, 
demonstrating the Brownian dynamics algorithm. Collisions with 
an obstacle are indicated by blue circles. Arrows indicate the di- 
rection and magnitude of the velocity after a collision (blue/dark) 
and after drawing a new velocity from the Boltzmann distribution 
(green/gray). 



the solvent's friction coefficient and temperature yield the co- 
efficient Dq of free diffusion. It may be combined with an in- 
trinsic length scale in the system, e. g., the obstacle radius a, 
yielding a microscopic time scale to :— a 2 /Dq. All averages 
for Brownian particles will refer to the canonical ensemble 
throughout the article. 

The overdamped dynamics of the tracer may be described 
by an effective stochastic force, caused by incessant colli- 
sions with solvent molecules. We have extended the above 
simulation algorithm by taking advantage of a coarse-grained 
scheme, where the ballistic trajectory is frequently interrupted 
by an equilibration with the solvent. After each fixed time 
interval tb a new velocity is assigned to the tracer particle, 
randomly drawn from the normalized Boltzmann distribution, 



p v (v) oc exp 



2fc B T 



(25) 



Any hydrodynamic interactions are neglected. Such a sim- 
ulation scheme for Brownian particles with hard-core inter- 
actions was carefully tested recently ll62Tl : a similar approach 
was already studied by Alley HU], where the equilibration was 
restricted to the instants of particle collisions. Fig.[3]exempli- 
fies the motion in two dimensions, demonstrating the velocity 
changes along the trajectory caused by events of collision and 
equilibration. 

Without obstacles, the algorithm yields a tb -dependent dif- 
fusion coefficient, cf. Eq. @, 



(v 2 ) 



TB 



2d 



k B T 
2m 



tb- 



(26) 



The obtained motion is obviously only diffusive for times 
larger than tb. To model Brownian motion at all physically 
relevant time scales, one has to choose r B smaller than the 
shortest time scale of the system; in the present problem, this 
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FIG. 4: (Color online) Convergence of the macroscopic diffusion 
coefficient of the Brownian particle for decreasing time interval tb at 
different densities. The limits D(tb — ► 0) are estimated by averaging 
over results for tb < 0.15 v a. 



is the average time interval r c between collisions with ob- 
stacles. The dependence of the diffusion coefficient on tb 
is shown in Fig. [4] We consider the value of tb sufficiently 
small when the macroscopic diffusion coefficient does not de- 
pend on tb anymore, resulting in tb « 0.3r c at criticality, 
n* = 0.839. 

Considering the velocity autocorrelation function ip(t) = 
(v(0) • v(f)) jv 2 for the combined algorithm of Brownian 
dynamics and ballistic collisions, one would naively expect 
that all correlations in the velocity vanish immediately for 
times t > tb. This seems reasonable, since a new velocity 
drawn at random should be uncorrelated to the previous value. 
The presence of excluded volume, however, induces correla- 
tions for particles sufficiently close to an obstacle if the Brow- 
nian update interval and the mean interval between collisions 
are comparable, tb r* r c . As a consequence, ip(t) shows pe- 
riodic anticorrelations for times t > tb , with the periodicity 
interval tb and an exponentially decaying amplitude. This ef- 
fect is demonstrated in Fig.[5]fbr several different values of tb. 
Since the focus of this work is on the long-time properties, we 
defer a more detailed explanation to Appendix lAl 



V. SIMULATION RESULTS 

We have simulated trajectories of ballistic and Brownian 
tracer particles over a wide range of obstacle densities, above 
and below the localization transition (Fig. |6). The mean- 
square displacement exhibits three distinct regimes: at short 
times, t -C to, transport is not hindered by the obstacles, and 
the intrinsic dynamics of the particle is observed, either bal- 
listic or Brownian. In both cases, this regime is followed by a 
regime of anomalous transport to be discussed in detail later. 
At large time scales, t ^> t x , and below the localization tran- 
sition, generic diffusive behavior is recovered. Above the crit- 
ical density, the particles are trapped and 8r 2 (t) is bounded 
by I 2 , the mean-square cluster size. At the critical density, the 
dynamics becomes neither diffusive nor trapped, and transport 
remains anomalous for all times. 




0.01 0.1 % c 1 10 



Time t/v a 

FIG. 5: (Color online) Velocity-autocorrelation function tp(t) ob- 
tained from the Brownian dynamics algorithm, at density n* = 0.6. 
Simulation results are shown for four different time intervals tb and 
for the limiting case tb — > oo, equivalent to a ballistic particle. ip(t) 
does not vanish immediately for t > tb, but oscillates with frequency 
1/tb. On the a>axis, the mean collision rate r" 1 is also indicated, 
varying inappreciably with tb. Inset: magnification of the velocity 
anti-correlations; vertical gray lines indicate tb of the data with the 
corresponding dashing pattern. Note that the chosen tb are multiples 
of the smallest one. 



A. Ballistic particles 

Let us discuss ballistic particles first. Diffusion coeffi- 
cients were extracted from the mean-square displacements, 
D = dtSr 2 (t — ► oo)/6. They vanish as a critical den- 
sity n* is approached from below according to a power law, 
D ~ \n* — n*]^. The exponent fi could be fitted to these 
data; such a fit, however, depends on a simultaneous determi- 
nation of the critical density n*. We will rather use a the- 
oretical prediction for fj, that relies on the mapping to ran- 
dom resistor networks; therewith, a first test of the validity 
of the mapping is obtained. Eq. (fl6l i tog ether with the re- 
sult a = 1/2 from Machta and Moore Eq. ®, provides 
\x = v + 2 w 2.88. This prediction is clearly corroborated 
by our results: in Fig. [7] the diffusion coefficient obeys the 
power law over five orders of magnitude, with a deviation of 
less than 15%. 

The knowledge of yu allows to fit the critical density quite 
precisely by means of a rectification plot, showing D 1 ^ 
against n* on a linear scale; see Fig. [7J:. We obtain n* = 
0.839(4), which coincides with the percolation threshold of 
the void space, n ; e[c = 0.8363(24) H gl B . The in- 
terpretation is that the particle will eventually squeeze through 
any gap, no matter how narrow. There are no regions on a 
cluster which are too improbable to be visited by the particle 
after an infinitely long time. This means furthermore that the 
particle will diffuse as long as its surrounding phase space is 
connected with infinity. 

In the localized regime, n* > n*, the long-time limit of the 
mean-square displacement directly yields the mean cluster ra- 
dius, Sr 2 (t) ~ £ 2 for t 3> t x . Then, the correlation length £ 
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FIG. 6: (Color online) Mean-square displacement Sr 2 (t) of the Lorentz model with (a) ballistic and (b) Brownian particles. The obstacle 
density n* increases from top to bottom; thick black lines indicate the long-time asymptote for anomalous transport at criticality, Sr 2 (t) ~ t 2 ^ z 
with z = 6.25. 




FIG. 7: (Color online) Critical power law behavior close to the localization transition; black closed symbols refer to ballistic particles, blue 
open squares to Brownian ones, (a) As the localization transition is approached, the diffusion coefficient D vanishes with exponent /i = 2.88. 
Units are va for ballistic and Dq for Brownian particles, (b) The localization length £ diverges with exponent v — /3/2 and can clearly be 
distinguished from the correlation length, £ ~ | e | ~~ " . (c ) Rectification of the diffusion coefficient close to the critical density. Fitting a straight 
line to the closest five data points yields the critical density n*. (d) Comparing D with its asymptotic behavior -D c |e| M on a semi-logarithmic 
scale reveals that the deviation from the power law is less than 10% for ballistic particles, except for the smallest data point. 



is easily accessible upon observing that the mean quartic dis- 
placement Sr 4 (t) := (AR(f) 4 ) scales for long times as l25ll 

5r 4 (t)~£ 2 £ 2 for i>i x ; (27) 

this relation is taken as definition for £. The extracted length 
scales I and £ diverge at the localization transition with expo- 
nents v — j3/2 = 0.68 and v = 0.88, respectively, according 
to Eqs. < f20b and ( fl~8l >; see Fig. |7J). The values for the expo- 
nents are taken from lattice percolation 15711 . and thus our data 
corroborate that the geometric exponents from lattices apply 
equally to continuum percolation 14611 . Further, the localiza- 
tion length I is identified with the mean-cluster size, as al- 
ready anticipated by our notation, and clearly contrasted with 



the correlation length £. Both length scales are relevant for 
the dynamics as proposed in Ref. 12511 . In particular, our re- 
sults corroborate the interpretation of the dynamic properties 
of the Lorentz model in terms of random walks on percolation 
clusters. 

There is a competing prediction /i = v + 3/2 ll42ll . which 
is also found in text books ll53ll . A third, maybe naive guess 
suggests ^i lat as 2.0, the universal value for lattice percolation 
in d = 3 l57ll . We will allow for these alternative predictions 
by a sensitive and unbiased test based on Eq. d22l , which re- 
lates the different values of /1 to different exponents z for the 
anomalous transport. The latter can be inferred directly at the 
critical density n* = 0.839 from the subdiffusive behavior of 
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FIG. 8: (Color online) Critical behavior of the mean-square displace- 
ment at n* = 0.839 ~ n c for ballistic and Brownian particles. 
Straight lines represent power laws t 2 ^ z , where z is related to [i by 
Eq. d22t . Three predicted values of /i are compared: the scaling rela- 
tions n = v + 2 l40ll and /i = v + 3/2 from Eq. l |16l > as well as 
the value on lattices, /i lat = 2.0. Our data are compatible only with 
the prediction by Machta and Moore l4fj|l . Inset: as a most sensitive 
test, the same data are rectified employing /i = v + 2; straight lines 
indicate the estimated long-time limits. The convergence for long 
times provides compelling support for this value of fj,. 



the mean-square displacement, growing as Sr 2 (t) ~ t 2 ^ z for 
long times, see Fig. [8] In the double-logarithmic plot, our data 
exhibit a slope manifestly smaller than expected from the al- 
ternative values for z. The data, however, collapse very well 
with a slope 2/z = 2/6.25 corresponding to /i = 2.88. For 
t > 10 5 er/u, we find a deviation of the mean-square displace- 
ment from this asymptotic behavior by less than 7% over time 
scales spanning 4 decades, see inset. In conclusion, only the 
value [i = 2.88 is consistent with our data, the other two can- 
didates can clearly be ruled out. Hence, the hyperscaling re- 
lation, Eq. ( [16) , with the value a = 1/2 holds for the Lorentz 
model. 



B. Brownian particles 

For Brownian particles, the different microscopic dynamics 
becomes immediately apparent in the short time regime of the 
mean-square displacement, see Fig. [(J. While ballistic parti- 
cles show a pronounced bending to diffusive or subdiffusive 
motion at the time scale of the collisions r c , Brownian parti- 
cles are hardly effected by the obstacles at low densities. At 
higher densities, the excluded volume induces a transient sub- 
diffusive regime again; the crossover from microscopic diffu- 
sion to the subdiffusive regime, however, is considerably more 
spread out. 

From universality arguments one expects that the diffu- 
sion coefficient vanishes again at some critical density with a 
power law that has the same exponent fi as in the case of bal- 
listic particles. Plotting D 1 ^ against density n* indeed yields 
a straight line, see Fig. [7J:. In addition, one infers that the crit- 
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FIG. 9: (Color online) The non-Gaussian parameter for Brownian 
particles yields finite values in the long-time limit, and it increases 
substantially as the localization transition is approached. Inset: at 
criticality, 02(f) grows with a power law in time; the predicted ex- 
ponent 0.097 is indicated by a straight line. 



ical density for Brownian and ballistic particles is the same, 
emphasizing that the transition is due to geometric rather than 
dynamic effects. Yet the asymptotic power law is markedly 
slower approached as in the ballistic case. 

The slow convergence for overdamped microdynamics is 
mostly pronounced in the mean-square displacement at criti- 
cality, see Fig.|8] For t w 10%, it is still more than 10% off — 
about a factor 100 slower than for ballistic particles. Never- 
theless, anomalous diffusion with the same dynamic exponent 
z is ultimately observed. 

The spatially heterogeneous character of the medium ren- 
ders the motion non-Gaussian, i. e., the distribution of particle 
displacements after a given time lag deviates from a Gaussian 
distribution. In supercooled liquids, dynamic heterogeneities 
have been quantified in terms of the non-Gaussian parame- 
ter JH, defined as S 



a 2 {t) 



3 5r 4 (t) 



(28) 



At moderate densities, the presence of the obstacles lets 012 (t) 
rise to values around 1 at intermediate times, until it decays to 
a finite value close to zero, see Fig. [9] This finite long-time 
limit is due to particles trapped in finite clusters; it diverges 
as criticality is approached. At criticality, the non-Gaussian 
parameter is predicted to grow with a power law, «2(i) ~ 
as a consequence of the competition between 
the localization length I and the correlation length £ 1I25I1 . Our 
data provide evidence for a strong increase of ct2(t) close to 
n*, and a double-logarithmic plot yields straight lines. The 
predicted exponent is very small, (3 /{2v — (3 + /i) = 0.097, 
and difficult to observe, but compatible with our data. 
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FIG. 10: (Color online) Scaling plots of the mean-square displacements for Brownian particles. Anomalous transport in the critical regime 
corresponds to a constant for t — > 0. Time is rescaled with the measured diffusion coefficients, (a) Scaling at leading order: as n* — » n*, 
data converge to 1 for small t, and an asymptotic data collapse may be anticipated. The amplitude A was determined from the inset of 
Fig. HI (b) Inclusion of corrections to scaling at leading-order, Eq. d31t with A_(f) = C = — 9.5^, yields excellent data collapse onto 
the scaling function 8r 2 _(i). Rescaled and corrected mean-square displacements for ballistic particles (C = — 0.8tg) are added in gray from 
Fig. 2b of Ref. 12^1 . The perfect match of Brownian and ballistic scaling functions substantiates universality of both systems. Inset: the 
corrections at criticality, n* = 0.839, decay with a power law for long times, see Eq. d32b . The exponent y = 0.34 is nicely corroborated over 
two decades in time, and the amplitude can be limited to the range —9 < Ct~^ v < — 15 (straight lines). 



C. Dynamic scaling 

A more stringent test of the universality of ballistic and 
Brownian particles compares dynamic scaling functions. The 
dynamic scaling ansatz for the mean-square displacement 
reads 



6r 2 (t;s)~t 2 / z 6rl(t/t x ) 



(29) 



for e — > 0, t 3> To, and with the crossover time scaling as 
t x ~ l z . Note that the relevant length scale for the mean- 
square displacement is the mean cluster size I rather than the 
correlation length £, due to the infinite life time of the perco- 
lation clusters. The subscript ± at the scaling functions refers 
to the sign of e, discriminating the different behaviors in the 
long-time limit. In the following, we will restrict the discus- 
sion to the diffusive regime, e < 0. The scaling function 
8r 2 _{t) with t := t/t x interpolates between anomalous dif- 
fusion at criticality and normal diffusion for long times; thus 



5r 2 _{t) 



A 

A rp-2/z 



for 
for 



0. 

oo. 



(30) 



The crossover time may be defined quantitatively by matching 
both regimes, At x z = &Dt x , implying that t x ~ £P/( 2_Z ) 
for e — > 0. For ballistic particles, the diffusion coefficient 
observes nicely the asymptotic law D ~ |e| M , see Fig. [7jJ, 
and plotting 5r 2 (t)/At 2/z vs. t\e\ 2u ~ l3+ ' 1 yields a satisfac- 
tory data collapse, see Fig. 2b in Ref. l25ll . For Brownian par- 
ticles however, our data for the diffusion coefficients are not 
yet in the asymptotic regime and show significant deviations 



from this behavior. Rescaling the mean-square displacements 
with t x in terms of e is not expected to lead to data collapse. 
Instead, let us rescale time with the measured diffusion coeffi- 
cients, t x oc ZJ-^-^+A 1 )/^ allowing for the deviations of D 
from its asymptotic behavior, see Fig.fTOb; the exponent eval- 
uates to (2v — (j + fj)/fj, ~ 1.47. All curves collapse in the 
diffusive regime (t 3> 1) by construction; but in the critical 
regime (i <C 1), the data fan out, and asymptotic convergence 
to a constant may only be anticipated. 

Such a behavior hints at corrections to the leading scaling 
behavior, which carry some fading reminiscence of the micro- 
scopic structure. One has to distinguish between analytic cor- 
rections, which depend, e.g., on the choice of the separation 
parameter, and universal corrections with non-integer powers 
(or logarithmic terms). In the case of ballistic particles l25ll . 
we have identified the leading correction as a universal power 
law, and excellent data collapse has been achieved in the crit- 
ical regime too by allowing for these terms. We have argued 
that the scaling ansatz for the mean-square displacement in- 
cluding the leading correction reads 



5r 2 {t;e) ~ t 2 / z 8r 2 ± {t) [l + t~ v A z 



(31) 



The exponent y is universal and will be discussed below. The 
correction function A±(i) is universal too, but unknown. At 
criticality, it reduces to a constant, A±(0) =: C, and the 
mean-square displacement obeys 



5r 2 {t;e = Q)~ At 2 '* (l + Cr v ) 



(32) 



The dynamic correction exponent y can be related to the 
static correction exponent £1 via the exponent relation 



yd* 



Qdf 



(33) 
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d w denotes the walk dimension introduced in Section [HI] and 
df = d — (3/v is the fractal space dimension. This relation 
has been derived by two of us within a cluster-resolved scal- 
ing theory for percolation l60ll . The main idea of the deriva- 
tion is to consider a propagator for tracer particles restricted 
to clusters of size s. Including an irrelevant parameter of the 
propagator, similar as in Ref. |25ll , generates the leading dy- 
namic correction, Eq. OTb . as well as the leading correction 
to the cluster size distribution n s , 



s- T [l + 0(s~ n )] for s -> oo, 



(34) 



with the Fisher exponent t = 1 + d/df. For lattice percolation 
in d = 3, the cluster size distribution has been analyzed with 
extensive Monte Carlo simulations, and the static correction 
exponent has been determined to = 0.64 ± 0.02 l67ll . As- 
suming that lattice and continuum percolation share the same 
geometric exponents, one calculates for the three-dimensional 
Lorentz model y = 0.34. 

For ballistic particles, the corrections are dominant in the 
critical regime, and they are well described by approximating 
A±(i) ps C. The corrections to the diffusion coefficient are 
encoded in A_ (t) for large t. They are not small for Brown- 
ian particles, and such a simple approximation will fail. But 
if these deviations are taken into account by hand as above, 
it seems reasonable to use A±(i) « C again. The value of 
C may be fitted at the critical density by means of Eq. d32l . 
see inset of Fig.fTOb. Our data corroborate the correction ex- 
ponent over two decades in time, but the determination of C 
is hindered due to statistical noise; the data permit a range 
9 < C < 15 (in units to = 1). Note that the correction ampli- 
tude is substantially larger than for ballistic particles, explain- 
ing the poor convergence to the long-time asymptote. Away 
from the critical density, we have rescaled the mean-square 
displacements again, taking into account the discussed cor- 
rections. Plotting 8r 2 (t)/At 2 t z (l + Ct~ y ) vs. time rescaled 
with the measured diffusion coefficients in Fig. [T0b . the data 
collapse perfectly onto the the scaling function 5r 2 _(t). 

The universality hypothesis predicts further that different 
systems are described by the same scaling functions, provided 
the latter encode the renormalization flow between the same 
two fixed points; in the present example, the fixed points are 
given by the critical point and the diffusive long-time limit. 
In Fig.fTOb, we have also included the collapsed mean-square 
displacements for ballistic particles from Ref. l25ll . An excel- 
lent match with the data for Brownian particles is observed, 
including the crossover regime over more than 5 non-trivial 
decades in time. Whence the amplitude A and the diffusion 
coefficients are determined (cf. Figs. Q and [SJ, no adjustable 
parameter enters the plot. One concludes that a single uni- 
versal scaling function 5r 2 _(t) describes the crossover from 
anomalous to normal diffusion for both ballistic and Brown- 
ian particles. 



D. Finite-size scaling 

The size of the simulation box L limits the size of the 
largest finite clusters; hence it acts as an upper cutoff on the 
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FIG. 1 1 : (Color online) The rescaled diffusion coefficients (for bal- 
listic particles) collapse on the scaling functions D±(x), which de- 
scribe the finite-size effects. For L <C £, the rescaled diffusion coef- 
ficient diverges asymptotically with exponent fi/v, see Eq. l !36b . 



correlation length. For very small system, L <C £, the cor- 
relation length becomes irrelevant, and the intrinsic (macro- 
scopic) length scale is given by the box size. For very large 
systems, L £, the finiteness of the box can be neglected. 
Observables and scaling functions are thus decorated with an 
additional parameter, L and L/£, respectively. We extend the 
scaling ansatz for the mean-square displacement, Eq. ( |29] l, to 



6r 2 {t; e, L) = t 2 ^Sr 2 {t£- z , L/£), 
and infer for the diffusion coefficient, taking t — > oo, 

D{e,L)=r^b ± {L/i) 
' L-^/ y for 



(35) 



O (e- L /«) 



for 
for 



and 
and 



(36) 



Note that in small systems, the diffusion coefficient becomes 
independent of e close to the critical density. In particular, 
diffusion is not blocked for n* > n* as long as the box size 
is smaller than the correlation length. Diffusing particles exist 
above n* c since the periodic boundary conditions turn large 
finite clusters that would exceed the box into infinite clusters. 
The predicted finite-size scaling of the diffusion coefficient is 
nicely observed in Fig. [TT] where data for different box sizes 
and densities above and below the threshold collapse into the 
two master curves D±. It follows that finite-size effects in the 
Lorentz model may be neglected as soon as L/( > 5, and 
most importantly, our foregoing analysis is not spoiled by the 
finiteness of the simulation box. 

The finite-size effects pretend a shift of the percolation 
threshold towards the localized regime l57ll . We have deter- 
mined the critical density from the axis intercept of the linear 
extrapolation of Z?(n*) 1 / A1 vs. n*, see Fig. [7]; and Fig.[T2lfor 
different box sizes. Rewriting Eq. (l36l l as 



D(e, L) 



-l/u 



D(eL 



(37) 
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FIG. 12: (Color online) Rectification plot of the diffusion coefficients 
of ballistic particles for various box sizes. Straight lines are fits to the 
data for n* < 0.82, their zeros yield the effective critical densities 
n*(L). Inset: extrapolation of n*(L) to infinite systems. 

with an analytic scaling function D, one infers that D^*) 1 ^ 
is shifted with increasing box size and that the effective criti- 
cal density n*(L) deviates from the true one as 

n*{L)-n* ~L- 1/v . (38) 

Plotting the effective critical densites n*(L) against L~ l l v in 
the inset of Fig.[T2]indeed results in a straight line. The extrap- 
olation to infinite systems allows a more precise determination 
of the critical density, we obtain n* = 0.837±0.002. We have 
checked that the result is robust against a possible error in the 
critical exponent v of the order of 0.01. Therewith, our pre- 
vious estimate from Section [V] based on a fixed box size, is 
improved and agrees very well with the result from Rintoul 
& n* = 0.8363 ±0.0024. 

VI. CONCLUSION 

The dynamics of ballistic and Brownian particles was in- 
vestigated in a heterogeneous environment close to the local- 
ization transition. It has been demonstrated that both systems 
share the same phenomenology on scales where the micro- 
scopic details are not resolved anymore. We have further 
corroborated our previous findings that the localization tran- 
sition is induced solely by a change in the topology of the 
medium: the percolation transition of the void space. Most 
importantly, our data substantiate that both Brownian and bal- 
listic particles in a percolating medium belong to the same 
dynamic universality class. Taking into account the leading 
corrections to scaling, we have extracted the universal scal- 
ing function for the crossover from anomalous, subdiffusive 
transport to normal diffusion. The crossover is found to vary 
remarkably slowly, spanning at least 5 decades in time. In 
addition, it is an interesting observation that the asymptotic 
regime is much slower approached by Brownian particles; this 
statement refers equally well to the asymptotic behavior of the 



mean-square displacement and to the suppression of the diffu- 
sion coefficient upon approaching the localization transition. 
As a consequence, to observe the genuine asymptotic power 
law of the anomalous transport is an experimental challenge; 
very large time windows are required to distinguish universal 
behavior from transient crossover phenomena. 

Although the divergent length scales I and £ cover only a 
decade in the investigated parameter regimes, we find anoma- 
lous transport already over many decades in time. These find- 
ings suggest that in experimentally accessible time windows, 
it is likely to observe apparent power laws with varying ex- 
ponents, even if the spatial heterogeneities extend just over a 
decade in length scale. In particular, the Lorentz model pro- 
vides a generic mechanism for anomalous transport in spa- 
tially heterogeneous media. 

As mentioned in the introduction, protein transport in cells 
is anomalous due to the macromolecular crowding. For pos- 
sible applications to such highly complex systems, the robust- 
ness of the presented scenario has to be discussed. Since the 
constituents of a cell interact certainly not via a hard potential, 
it is tempting to replace the obstacles by soft spheres. As long 
as the kinetic energy of the (ballistic) particle is fixed, the dy- 
namics of a single hard particle, however, can be mapped one- 
to-one to a soft potential. A canonical ensemble seems more 
natural for a soft potential, but subsequent averaging with the 
Boltzmann weight smears out the transition and the critical 
properties since some fast particles can always overcome the 
obstacles. We have checked for the two-dimensional Lorentz 
model that this approach indeed reproduces simulation data 
for the diffusion coefficient 16811 . 
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APPENDIX A: SHORT-TIME VELOCITY CORRELATIONS 

In Section IIV Bl we have pointed on the effect of non- 
vanishing velocity correlations for t > tb in the Brownian 
dynamics simulations. Here we explain this behavior quanti- 
tatively within a toy model. 

The essential components for the explanation are, first, a 
random walk which is continuous in space but evaluated on a 
fixed time grid, and second, a hard wall restricting the avail- 
able volume through ballistic reflections. In the most simple 
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configuration, we have a point-like particle performing a two- 
dimensional random walk between two parallel walls. 

The calculation of the velocity autocorrelation function 
ip(t) is elementary for the first two time steps tb. Choos- 
ing a coordinate system with the walls parallel to the y-axis, 
we parametrize the velocity during time step n by an angle 
(p n to the x-axis, v„(t) = v„(tt,tb < t < (n + 1)tb) = 
v(cos (p n , sin if n ). There are four cases to be distinguished, 
depending on whether the particle collides with the wall or 
not during the first and the second time step, respectively. 
As an example, when we have collisions in both cases and a 
fixed absolute value v of the velocity, the correlation between 
the initial velocity vo(0) and the velocity one time step later, 
vi(£ > tb), reads 



where L is the distance between the walls, and the boundaries 
of the inner integral are defined by 



COS ifo 



-y; cos^o 




At > 1. 

T B — 



The same applies to the other three cases to be considered. 

The results obtained from this calculations are in excellent 
agreement with simulations of the toy model, see Fig. [13] In 
principle, %p(t) can be obtained for longer times t > 2tb too, 
but the calculations are quite involved. 



(vo(0) • Vi(tb + At)) = v (simpo sin tpi — cos ipo cos <pi). 

(Al) 

By geometric reasoning, one can reduce the averages to a dou- 
ble integral to be solved numerically, 



(v (0) • vi(r B + At)) 

l 4>o 

2v 3 t s 



Ltt 2 



dy jd(p' cos Poyl - (y) 



ipo 



y) (cos^ + y) 2 , 



(A2) 
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FIG. 13: (Color online) Velocity-autocorrelation function for the 
two-dimensional motion of a point-like particle moving between two 
walls oriented parallel to the y-axis. The walls are separated by 
a distance L in s-direction. The particle is reflected ballistically 
at the walls, and the direction of the velocity is changed randomly 
in intervals tb. The inset displays a magnification of the regime 
tb < t < 2tb on the linear scale. The analytic curve responding 
to Eq. &A2\ is shown for L = 2vtb ; it matches perfectly the simula- 
tion results. 
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